Conformal properties of ID quantum systems with long-range interactions 
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We investigate conformal properties of a one-dimensional quantum system with a long-range, 
Coulomb-like potential of the form with a > 0. We compute the conformal anomaly c as 

function of a. We obtain c = for a < 1 (Wigner crystal) and c — 1 for a > 1 (Luttinger 
liquid). By studying the scaling regime of the system when it is deformed by the inclusion of a 
term t p(x), (with p{x) a scaling operator), we show that, for cr > 1, the correlation length ^ gets an 
additional interaction dependent factor. This leads to a necessary redefinition of ^ in order to avoid 
an unphysical dependence of the central charge on the coupling constant. 
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I. INTRODUCTION 

The principle of conformal invariance (CI) at a criti- 
cal point is one of the most powerful tools for our un- 
derstanding of one-dimensional (ID) quantum systems^. 
This invariance constrains the possible behavior of a 
given critical theory, allowing a characterization of uni- 
versality classes in terms of a dimensionless parameter c 
(the conformal anomaly) which is the central charge of 
the quantum extension of the (infinite-dimensional) con- 
formal algebra. The central charge gives a measure of the 
sensitivity of a critical theory when one changes the space 
geometry. It has also been related to the low-temperature 
specific heat in quantum chains^. As a striking achieve- 
ment, in the context of classical two-dimensional statis- 
tical mechanics models, conformal invariance led to the 
derivation of all critical exponents corresponding to the- 
ories with c < l'^-. When the models posses an additional 
continuous symmetry a classification of critical models is 
also possible for c > 1—. 

CI is obtained by assuming translational, rotational 
and scale invariance, together with short-range interac- 
tions. If any of these conditions are not verified its valid- 
ity can not be assured. However, the concepts of CI have 
been extended in several directions. Cardy, has derived a 
relationship between the second moment of the correla- 
tion function of a scaling operator away from criticality 
(in the scaling region) and the central charge and the 
scaling dimension corresponding to the critical theory^. 
CI has been also applied, as a principle, to ID quantum 
systems with inverse-square long-range interactions^ i. 
In particular, the continuum limit of the 5=1/2 Heisen- 
berg chain with inverse-square exchange was associated 
to a c = 1 SU{2) Kac-Moody theory^ and Sutherland's 
model with interaction g/r'^ ^ was analyzed through fi- 
nite size scaling^ . Finite size corrections to the ground 
state led to a g-dependent central charge, whereas the low 
temperature expansion of the free energy yielded c — 1, 
which is consistent with the scaling dimensions obtained 
for density operators and momentum distributions. The 



reason for this discrepancy has not been fully explained, 
although it was suggested that it might be related to the 
finite-size approximation to the long-range potential^ and 
to the possible differences between time and space corre- 
lations in long-ranged modelaiS. 

On the other hand, in the context of low-dimensional 
condensed matter systems, the interest in long-range in- 
teractions is not a purely academic one. For instance, 
in quasi-one-dimensional systems one usually considers 
contact electron-electron interactions assuming that the 
presence of adjacent chains effectively screens the ba- 
sic Coulomb interaction. However, as the dimensionality 
of the system decreases charge screening effects become 
less important and consequently the long-range interac- 
tion plays an observable roleii. The effect of power- 
law long-range interactions is also of interest for the 
study of antiferromagnetic spin-chains^^i, in which pos- 
sible experimental realizations might be found, for in- 
stance, in systems with Ruderman-Kittcl-Kasuya-Yosida 
interactionsi^. 

The purpose of this paper is to explore the interplay 
between the principle of conformal invariance and the ef- 
fect of long-range, distance-dependent interactions. By 
performing analytical computations of the central charge 
c in a continuous theory, we shed light on the effect of 
anisotropy on the conformal properties of a ID quan- 
tum system. In particular, in connection to the above 
mentioned discrepancy, we will explicitly show how the 
asymmetry between long-time and long-space behaviors 
of the correlations must be treated in order to obtain 
sensible results for the conformal anomaly. 

Most of the previous calculations of c for long-ranged 
systems have been done on the basis of a discrete model 
and for specific values of a. The exception is the impor- 
tant paper by Mironov and Zabrodir.^'^, who were able 
to determine scaling dimensions and correlation functions 
for a general class of potentials. However, in their formu- 
lation it was not possible to specify the conditions to be 
satisfied by the potentials in order to assure the absence 
of a gap. Moreover, the value of c was not computed but 
inferred from the form of some scaling dimensions. 
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II. THE MODEL AND THE LONG-RANGE 
POTENTIAL 

We will study a bosonic non local (l+l)-dimensional 
quantum field theory with Lagrangian density 

C^=^-{d^ct^{x)f + ]^ J d^yd,c^{x)V{x~y)d^c^{y). (1) 

This model is the bosonized version of a Thirring- 
Tomonaga-Luttinger fermionic model containing only 
forward-scattering interactions-'^^. It describes charge- 
density oscillations with sound velocity v = -^/l -I- V{pi), 
where V{pi) is the Fourier transform of the distance- 
dependent potential. Although the model under consid- 
eration is not Lorentz-covariant, we shall use, for con- 
venience, a notation which is reminiscent of relativistic 
problems {cPx = dxodxi, d'^p = dpodpi). 

We will consider an instantaneous, power-law decaying 
potential of the form 

V{x)=g5{xo)ixj+a')^, (2) 

where a is an arbitrary positive real number and a is a 
regulator parameter that can be associated to the lattice 
spacing. 

Since we shall deal with all positive values of a we will 
need the Fourier transformed potential which depends on 
a as follows: 

(3) 

for a > I. Here -ft'(i-(j)/2(a | Pi |) is a modified Bessel 
function. Note that, in this case, the potential is renor- 
malized. This renormalization is associated to the diver- 
gence of the integrand in the Fourier integral, for xi — 0. 
This divergence is not present when a < 1. Defining the 
renormalized coupling constant gjj = {a)^~'^g and tak- 
ing into account the behaviour of the Bessel function for 
small arguments we get 

V{pi)^f>{9R,a), (4) 

with/>(gH,a)^ ^-v^ff;-^)/^) . 

For CT < 1 the Fourier-transformed potential reads 

V{pi) = f<{9,a)\p,\<^^-'\ (5) 

with/<(g,.)^ .^v^%7J,y--)/^\ 

As we shall see, the Coulomb potential correspond- 
ing to (T = 1 is a special case that can be considered 
as the frontier between two regions of qualitatively dif- 
ferent conformal behaviors. This case was investigated, 
from another point of view, by Schulz, in his work on 
Wigner crystal-like formation in a Luttinger systemic 
(See also ref. 17 for an explicit computation of the elec- 
tronic single-particle function). For a = 1 one has 

Vipi) = 2gKo{a\pi\) « -2 g Ha\pi\). (6) 



III. COMPUTATION OF THE CONFORMAL 
ANOMALY 

Exactly as it was done in Refs. 9 and 10 with the 
Sutherland's model, for the special case a = 2, in this 
paper we shall employ the model given by Q (for arbi- 
trary real and positive a) as a testing ground to exam- 
ine the consequences of imposing CI in the presence of 
anisotropic interactions. Our idea is to compute the cen- 
tral charge following two different routes. Firstly, we will 
compute the low-temperature behaviour of the specific 
heat, from which we can read the value of c, according to 
the results of Refs. 2. Secondly we will perturb the crit- 
ical Lagrangian by adding a simple scaling operator 
p{x) that takes the system away from criticality. Then, 
by considering the second moment of the two-point cor- 
relations of p{x) in the scaling regime, we can obtain the 
value of the conformal anomaly by using the result de- 
rived by Cardy^: 

/d^.l.HK.Mo)).^^-^,^!-^, (7) 

where is the scaling dimension of p and t is the cou- 
pling constant of the added term that spoils CI {t cx 
cx (T — Tc)). Although somewhat indirect, this 
method is interesting because it involves both large spa- 
tial and temporal dependence of the correlations. This 
fact will enable us to unravel some subtleties involved 
in the application of CI concepts to anisotropic models. 
Let us point out that, eventhough the above equation is 
valid for any scaling operator, p{x) is the density opera- 
tor throughout this work. 



A. First method: free-energy calculation 

In order to compute the low temperature behavior 
of the specific heat we study iQJ in the imaginary-time 
formalisml^. As usual, we build the partition function 
Z as a functional integral extended over the paths with 
periodicity conditions in the Euclidean time variable xq- 
4>{xq +(3^x1) = (f>{xo,xi), where /3 = y,^:^ = 1- After 
some standard algebraic steps (see Ref. 15 for a detailed 
similar computation in a fermionic model) the Helmholtz 
free energy !F — — In Z can be expressed as 

^ = + i r ^ln(l-e-^l''^l''(^'^)), (8) 
(3 J_oo 27r 

where J^q is the zero-point, vacuum free energy. Now 
we have to evaluate the integral and read the value of c. 
This is easily done for cr > 1 , since v — y^l + f>{gR, a) — 
v{gii, a) is independent of the momentimr. The result is 

6p^w(5fl,cr) 
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which means that c = 1—. For a < 1 the computation is a 
httle bit more subtle because the sound velocity depends 
on pi. However, taking into account that one is inter- 
ested in the large-/? behaviour of !F, the integral can be 

approximated, yielding a result proportional to P 
For a < 1 the exponent > 2, from which one con- 
cludes that c = 0. Conformal field theories with c = 
have been previously found in the study of second order 
phase transitions in the presence of quenched disordenSii. 
Some of their interesting properties have been analyzed 
quite recentljiSi. 



and try to compute the left hand side of ^ . This compu- 
tation involves correlations slightly away from criticality, 
in the scaling regime in which both |a;i| and the correla- 
tion length ^ are very large, but with |a;i| << ^. Besides, 
both spatial and temporal dependence of the correlations 
are needed. In fact, regarding the possible application of 
CI to an anisotropic model such as the one considered 
here, this is a crucial point. Indeed, it is very instructive 
to compute again {p{x) p(0))o, not in the equal-time case 
but for Xq,Xi 0. This is quite easily done for a > 1. 
The result is 



B. Second method: deformation of a conformal 
field theory 

Let us now undertake the computation of c following 
the second method described above. We choose as per- 
turbing operator p(x) — di4>(x). In order to determine 
the scaling dimension we compute the critical two-point 
correlation. A straightforward manipulation yields 



{p{x) = 



pI+pI v'^ipi) 



(10) 



As a first step we evaluate the scaling dimensions by 
considering equal-time correlations. Performing a long- 
distance approximation in the above integral we get 



{p{x,) p{0))o ^ 
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(11) 



for cr > 1, from which we obtain Ap — 1. Note that in 

this case v = \/l + f>igR,cr) = v{gR, a). 

The analogous calculation for ct < 1 gives a lengthy ex- 
pression in terms of a confluent hypergeometric function. 
Analyzing the asymptotic long distance limit we obtain 



(p(a;i)p(0))o 



7r(5- 0-) y/f<{g,(j) \xi\'^ 



(12) 



from which we get Ap ~ 1 + Note that both results, 
for cr > 1 and a < 1 coincide for a = 1, which suggests 
that Ap is continuous for all a > 0. We have explicitly 
verified the correctness of this assertion: 



(p(x)p(0))o = ^^Er fed-- 2ln^], (13) 

where Erfc stands for the complementary error func- 
tion. In the long-distance limit this expression behaves 
as |a;i|^^, i.e. Ap ~ 1 for a — 1. We will use this result 
later. 

At this point we start considering the modified, non 
critical theory described by 



(14) 



{p{x)p{Q))o = 



-\xa\ v/\xi\ (■gbol v/\xi\ 



(15) 

Performing an expansion of the numerator around xq — 
we reobtain our previous result for the equal-time corre- 
lator, equation Hll|) . The important observation here is 
that although one can define a scaling dimension which is 
uniquely determined for large distances and long times, 
the precise functional form of the correlation depends on 
the direction in which large scales are observed. The in- 
teresting point is that in contrast to previous calcula- 
tions, our procedure is able to take into account this ex- 
pected anisotropy. Since the integrand in the left hand 
side of (|7|) involves the general (not equal-time) behavior 
of the correlator, our result for c will reflect the contri- 
butions coming from different directions of space-time. 
In view of these remarks we now compute the left hand 
side of I?!). In our calculation the correlation length is 
associated to spatial fluctuations only. This parameter 
is used as an infrared cutoff for the spatial integral. For 
(7 < 1 we get 



d^x \x\^{p{x)p{0))t = Ci r + C2 + ^3^2", (16) 



where Ci, C2 and C3 are constants depending on g and 
a. Since t ex. one sees that c — > 0, in agreement with 
the result obtained by computing the free energy ©. It 
is more illuminating the case a > 1 which gives 



(17) 



Since we know from ® that c = 1, we obtain an expres- 
sion for the correlation length in terms of t and v. The 
extra dependence of ^ onv is due to the anisotropy of the 
system, although the precise form of this dependence is 
a consequence of our regularization prescription. Please 
note that if we assume that ^ is a function of t only, 
the above equation predicts a conformal anomaly that 
depends on the long-range potential through the sound 
velocity, in analogy with the results previously obtained 
in Refs. 9 and 10 (for a different model and for a — 2). 
This situation reflects the fact that the correlation length 
is defined up to a multiplicative non- universal constant. 

For the particular value a = 1 the above procedures 
are not easy to implement due to the peculiar form of 



4 



the potential ©. However, based on the results and 
it is straightforward to conclude that c = 1 for u = 1 

also. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have presented the first explicit and 
analytical computation of the conformal anomaly in a 
system with long-range Coulomb-like interactions. For il- 
lustrative purposes we have considered the simple bosonic 
QFT defined in However, this model is inter- 

esting by itself since it is a bosonized version of the 
Tomonaga-Luttinger model with a distance-dependent 
forward-scattering potential. We obtained the central 
charge c as function of a. We predict c = for a very 
weakly decaying potential corresponding to cr < 1. Ac- 
cording to the results of Schulz^^ this corresponds to the 
regime in which the system behaves like a ID version of 
a Wigner crystal. For ct > 1 we obtain c = 1, in agree- 
ment with "common knowledge" on Luttinger liquids'-^. 
However, in previous computations the value of c was in- 
ferred from the scaling dimensions, whereas we are now 
providing a direct determination for a whole family of 



well specified Coulombian potentials. In contrast to all 
previous studies of CI in long-ranged models we have 
explicitly computed correlators for different spatial and 
temporal points^ which allowed us to capture the effect of 
anisotropy on c in a more complete fashion. Moreover, by 
considering as a mathematical tool a deformation of the 
critical theory given by H14|l and using a quantitative pre- 
diction relating c with properties of the deformed model 
in the scaling regime, we showed how, for cr > 1, the cor- 
relation length acquires a finite, interaction dependent 
renormalization factor. If one naively ignores this factor, 
an unphysical, interaction dependent conformal anomaly 
is obtained. 
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